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Origin of the optical second-harmonic generation in spherical gold nanoparticles: local 

surface and non-local bulk contributions 

G. BachelieiO and J. Butet, I. Russier-Antoine, C. Jonin, E. Benichou and P.-F. Brevet 
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The second-harmonic generation of 150 nm spherical gold nanoparticles is investigated both ex- 
perimentally and theoretically. We demonstrate that the interference effects between dipolar and 
octupolar plasmons can be used as a fingerprint to discriminate the local surface and non-local 
bulk contributions to the second-harmonic generation. By fitting the experimental data with the 
electric fields computed with finite-element method (FEM) simulations, the Rudnick and Stern pa- 
rameters weighting the relative nonlinear sources efficiencies are evaluated and the validity of the 
hydrodynamic model and the local density approximation approaches are discussed. 
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I. INTRODUCTION 



Nanoscience and nanotechnology have focused a wide 
interest in the past decades on noble-metal nanostruc- 
tures in view of the unique optical properties offered 
by the surface plasmon resonances (SPR). Their under- 
standing is nowadays reaching one of its ultimate devel- 
opment with the ability to combine numerical simula- 
tions, 3D transmission electronic microscopy, optical ex- 
tinction measurements or electron loss spectroscopy at 
the single particle level.—"— In addition, new opportuni- 
ties in plasmonics and meta materials have been recently 
opened with the coupling electric and magnetic dipoles 
in split-ring resonators or the tailoring of the SPR using 
Fano interferences. ■4'fi Though, one of the main emerging 
challenges lies in the dynamical, time-resolved control of 
the SPR. In this context, one appealing route is the use 
of nonlinear optics driven by femtosecond laser pulses. 
The understanding of the nonlinear processes at play in 
plasmonic nanostructures is therefore becoming a very 
active topic. 

The second-harmonic generation (SHG), whereby two 
photons at the fundamental frequency are converted into 
a single photon at the harmonic frequency, has quite a 
long history as it will soon celebrate its 50 years anniver- 
sary with the pioneer works of Franken et al.^ The ori- 
gin of the nonlinear sources has early been established 
as arising from the breakdown of the centrosymmetry at 
the metal surface (local response) or from field gradi- 
ents inside the bulk (non-local response). Two theoreti- 
cal approaches were pursued to account for their relative 
intensity: the analytic hydrodynamic model introduced 
by J.E. Sipe et al<i and the density functional approach 
proposed by A. Liebsch.® The later was found to be the 
most accurate description, although a complete theoreti- 
cal framework involving both intraband and interband 
transitions is still missing. Noble-metal nanoparticles 
have triggered a renewed interest in SHG. Especially, the 
size dependence of the nonlinear efficiency, the effects of 
resonance, composition and chirality of the nanostruc- 



tures have been addressed i2r2^ Single particle sensitiv- 
ity has also been reached in different configurations.—""— 
Analytical models have been developed for simple geome- 
tries such as spheres and cylinders.—""— If these impor- 
tant works mainly derived the selections rules and there- 
fore the excited and radiating multipolar surface plasmon 
modes, they did not quantitatively evaluate the relative 
contributions of surface and bulk sources to the SHG, 
which therefore remains a largely open question. Few 
recent works involving numerical simulations have sug- 
gested that the bulk currents are indeed a fundamen- 
tal contribution to the nonlinear response in noble-metal 
nanostructures ^ 

In this work, we investigate both experimentally and 
theoretically the nonlinear optical properties of 150 nm 
spherical gold nanoparticles in solution. We demon- 
strate for the first time that the interference effects be- 
tween dipolar and octupolar plasmons, which has been 
addressed recently22i, can be used as a fingerprint to dis- 
criminate the local surface and non-local bulk contribu- 
tions to the second-harmonic generation. In particular, 
the Rudnick and Stern parameters, weighting their rel- 
ative efficiencies ^21 are evaluated by fitting the experi- 
mental data with the electric fields computed with finite- 
element method (FEM) simulations. 



II. RESULTS AND DISCUSSION 

A. Local surface and non-local bulk contributions 

For isotropic and centrosymmetric materials, the 
second-order surface susceptibility tensor can be dras- 
tically reduced to only three independent components: 
X-L-L-Li XJ-IIII ^-nd X||||± where _L and || stem for perpen- 
dicular and parallel to the surface, respectively. However, 
from both theoretical and experimental point of vicw,^5i 
the XJ-IIII component only weakly contributes to the SHG 
of noble-metals. Hence, the two dominant surface non- 



linear polarizations can be written as 



P. 
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where the electric fields are evaluated inside the metal. 
These polarization vectors lead to surface currents lo- 
cated just outside the metal^ and given by jsurf — 
dPsurf/dt. From general considerations, the non-local 
bulk polarization may also have different contributions. 
However, in the special case of noble metals probed by 
a single fundamental field, it reduces to the following 
expressions^ 



bulk 



(r,2w) ^76„,fcV.[E(r,w).E(r,Lj)] 



(3) 



where jbuik is the bulk susceptibility. The associated 
bulk current is given by jbuik = dFbuik/dt. Hence, the 
knowledge oi x±±±, X\\\\± ^-^d ^buik allows in principle to 
completely determine the nonlinear optical properties of 
noblc-metals nanostructures. Following Sipe et alJ these 
parameters can be linked to the adimensional Rudnick 
and Stern parameters'-^ through the relations 
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(6) 



where e and m are the electron charge and mass and 
eoer(w) stems for the dielectric function of the metal at 
the frequency w.-SS In the framework of the hydrodynamic 
model, ^ the Rudnick and Stern parameters have very 
simple expressions: a = 1, 6 = — 1 and d = 1. However, 
if both tangential and bulk currents are well described 
by macroscopic parameters related to the dielectric func- 
tions, the normal current is driven by the local distribu- 
tion of the electronic density over the Fermi length below 
the metal surface. Therefore, one has to go beyond the 
hydrodynamic model in order to properly evaluate the 
parameter a. In particular, the latter was found to be 
larger than 1, frequency dependent and complex valued 
using local density approximation applied to conduction 
electrons only.- It is also important to note that both 
the hydrodynamic model and the local density approxi- 
mation approach were designed for planar surfaces and 
not for spherical particles. However, for particle radii by 
far larger than the Fermi length, the surface curvature ef- 
fect can be assumed to be negligible, as well as quantum 
size effects. 




Figure 1. Spatial distribution of the nonlinear currents cal- 
culated for the bulk contribution (a), the tangential (b) and 
normal (c) surface currents for a 150 nm gold nanoparticle 
excited at SOOnm. The near-field map of the corresponding 
harmonic field amplitudes are given in panels d-f. 



B. Finite-element method simulations 

The nonlinear currents corresponding to Eqs.[T]|n]were 
computed for a 150 nm gold nanoparticle assuming a = 1, 
fe = — 1 and d = 1 (see Fig. la-lc). The fundamental elec- 
tric field was obtained through finite-element simulations 
(FEM) using the scattered field formulation and perfectly 
matched layers in order to avoid any unwanted reflec- 
tion at the simulation boundaries.'^ In a second step, 
the near-field distribution of the harmonic field was cal- 
culated solving the Maxwell's equations within the weak 
formulation, allowing to incorporate both surface and 
bulk sources. As shown in Fig. Id-lf, the obtained near- 
field maps around the nanoparticle do not strongly de- 
pend on the nonlinear source, except in terms of overall 
intensity. The latter is directly driven by the parame- 
ters a, b and d but also by the evanescent and therefore 
localized character of the harmonic field. The knowl- 
edge of the far-field SHG intensity is a key ingredient to 
compare simulations with experimental measurements. 
It was computed from the near-field distribution of the 
harmonic field using the Stratton-Shu formulai^ 

Figures 2a and 2b show the SHG intensities simulated 
for an excitation beam (wavevector k) along the z axis 
and a collection at right angle (along the y axis) . The fun- 
damental electric field polarization angle with respect to 
the X axis is noted 7. For 7 = the polarization is along 
the X axis i.e. perpendicular to the scattering plane and 
7 = 90 deg corresponds to a polarization along the y axis, 
i.e. in the scattering plane. Depending on the selected 
polarization for the harmonic field, different patterns are 
observed. For a detection polarization perpendicular to 
the scattering plane, a four lobe pattern is obtained what- 
ever the origin of the nonlinear sources, arising either 
from the surface or from the bulk of the nanoparticle 
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Figure 2. Normalized far-field second-harmonic intensity 
computed using FEM simulations (a and b) and measured 
experimentally (c and d) as a function of the electric field 
polarization at the fundamental frequency. The polarization 
of the harmonic field is perpendicular (a and c) and parallel 
(b and d) to the scattering plane. The normal surface cur- 
rent (dashed-dotted lines), tangential surface current (dotted 
lines) and bulk currents (full lines) contributions are shown 
separately. 



(all curves are superposed in Fig. 2a). Invoking spherical 
harmonics characterized by the angular momenta {l,m) 
with —l<m<l, the latter property can be explained 
considering that i) the (k, E(r,aj)) plane and the plane 
normal to it are symmetry planes for the harmonic field 
allowing only even values for m, ii) the component of the 
field perpendicular to symmetry planes cancels out im- 
posing TO ^ for detection polarization perpendicular to 
the scattering plane and iii) the fundamental field con- 
tains only m = ±1 terms leading to polarization vectors 
in Eqs. [TJO having \m\ < 2. Hence, a quadrupolar pat- 
tern associated to to = ±2 is always observed, whatever 
the excited multipole given by /. This detection config- 
uration is therefore unlikely to discriminate the different 
contributions from the polarization measurements. 



C. Discrimination between the nonlinear sources 

Owing to the surface and bulk origin of the nonlinear 
currents, it might be tentative to use the size dependence 
of the SHG intensity to assess the relative contributions 
of the nonlinear sources. As a matter of fact, it was 
found experimentally that the SHG intensity scales as 
the squared surface for small gold nanoparticles whereas 
a squared volume dependence is observed for larger oneSfS 
However, unexpected dipolar patterns, normally forbid- 
den owing to the previous discussion, were revealed for 
small particles and for detection polarization perpendic- 
ular to the scattering plane. This was attributed to the 



centrosymmetry breaking of the particle shaped as later 
demonstrated with FEM simulations.— Indeed, all non- 
linear sources lead to the same size dependence of the 
scattered intensity, regardless of their surface or bulk ori- 
gin. This was verified with FEM simulations for spherical 
and deformed particles (data not shown here) , but can be 
directly understood from Eqs. [TEl the gradient appear- 
ing in the bulk polarization of centrosymmetric materials 
(Eq. [3]) introduces in spherical coordinates an additional 
term scaling as l/r and leading therefore to the same 
size dependence for the volume and surface polarization 
induced SH intensities. As a consequence, the size de- 
pendence cannot be used to separate bulk and surface 
contributions in centrosymmetric materials. 

Another strategy has therefore to be pursued in order 
to weight the nonlinear surface and bulk contributions in 
gold nanoparticles. This is provided by the interference 
effects between selected dipoles and octupoles we have 
recently demonstrated^^ for a detection polarization in 
the scattering plane. In this specific configuration, con- 
structive and destructive interferences are controlled by 
the incident polarization angle, leading to the varying 
intensities shown in Fig. [2]d. The key point here is that 
the interference pattern depends on the nonlinear source 
(surface or bulk). More precisely, if both surface terms 
have nearly the same polarization dependence for the 
SHG intensity, they clearly deviate from the bulk con- 
tribution for which the intensity weakly oscillates for a 
150 nm gold nanoparticle excited at 800 nm. 



D. Quantitative determination of the nonlinear 
sources efficiencies 



The experimental data shown in Fig. [5] were obtained 
using a modelocked Ti:sapphire laser tuned to a wave- 
length of 800 nm and delivering pulses of about 180 fs 
at a repetition rate of 76 MHz. The pulse energy mea- 
sured at the laser exit was 10 nJ. The laser beam was 
focused onto a quartz cell with a microscope objective 
(xl6, NA=0.32) leading to a beam waist at the funda- 
mental frequency of 5 ^m^ A low-pass filter was used 
in order to remove any residual light at the harmonic fre- 
quency generated prior to the cell. The SH photons are 
collected perpendicularly to the incident beam with a 25 
mm focal length lens (NA=0.5). The scattered photons 
at the fundamental frequency are removed by a high-pass 
filter placed before the monochromator. The polarization 
angle of the fundamental beam is selected with a rotating 
half- wave plate and the polarization of the SH photons is 
selected by an analyzer. The photon detection was per- 
formed by a sensitive cooled photomultiplier tube and 
the fundamental beam was chopped at 130 Hz allowing 
a gated photon counting regime in order to remove the 
background light. The colloid solution of gold nanopar- 
ticles dispersed in water (1.4 pM) were purchased from 
BBI International (average diameter of 150 nm with a 
standard deviation of i 



Comparing the calculated and measured signals (see 
Figs 2b and 2d), one can directly conclude that the nor- 
mal surface current, which is usually considered as dom- 
inating the nonlinear response, cannot account by itself 
for the observed pattern, in contrast with the results ob- 
tained in Ref. l23| for particle sizes smaller than 100 nm. 
Bulk and surface nonlinear sources being excited coher- 
ently, the corresponding SHG intensities shown in Fig. 2b 
cannot be summed up directly: the relative phases of the 
harmonic fields have to be taken into account. More pre- 
cisely, it has to be noted that these phases depend on the 
input polarization angle since different multipoles are in- 
volved. This key feature allows the discrimination of the 
two surface contributions, having otherwise nearly the 
same intensity pattern (see Fig. 2b). The experimental 
data were therefore fitted by 



IsHG =G'|(aEs„r/,±(r, 2w) + feE^^^/jKr, 2tj) 



(7) 



where a, b and d are the adjustable Rudnick and Stern 
parameters and G is a normalizing constant. Es„r/,j_, 



E 



surf, II 



and Ebuik are harmonic fields associated with 



the source terms given by Eqs. [T][3] as obtained from 
the FEM simulations in the far-field region. The inten- 
sities recorded for a harmonic polarization in and out of 
the scattering plane were simultaneously fitted using the 
same Rudnick and Stern parameters in order to fully ac- 
count for all experimental data. Owing to the normaliz- 
ing constant G, their absolute values cannot be extracted 
by the fitting procedure, but the relative amplitude and 
phase do. 



E. Comparison with existing models 

Different constraints were applied to the fitting Rud- 
nick and Stern parameters in order to check the valid- 
ity of existing models. We first considered the hydrody- 
namic modet where the Rudnick and Stern parameters 
satisfy a = d = —b = 1. Clearly, the obtained result 
cannot account for the experimental data, as shown in 
Fig. 3a. More advanced models were early proposed in- 
cluding a damping term to the hydrodynamic modeP^ or 
taking into account the spatial distribution of the elec- 
trons over the Fermi length below the surface within local 
density approximationJ^ All these models suggest that 
d ~ —b~ 1 and a is complex valued and exhibit res- 
onant behavior largely exceeding in amplitude the two 
other parameters. As a matter of fact, it was here again 
not possible to reproduce the experimental data by the 
fitting procedure (see Fig. 3b): additional minor peaks 
appear at 90 and 270 deg, in contrast with the experi- 
mental data. The fact that the damped hydrodynamic 
model and the local density approximation approaches 
can't accurately account for the SHG in gold nanoparti- 
cles is not surprising since they do not take into account 
the interband transitions that are resonantly excited at 
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Figure 3. Fits of the experimental data (only the harmonic 
polarization in the scattering plane is shown) using the com- 
plex valued harmonic fields, as computed with the FEM sim- 
ulations, and the following restrictions on the Rudnick and 
Stern parameters: (a) a = d = —b = 1, (b) a complex val- 
ued and d — —b = 1, (c) a complex valued, b real valued 
and d = 1 and (d) a complex valued, 6 = and d = 1. 
The scattered intensities corresponding to the normal surface 
current (dashed-dotted lines), tangential surface current (dot- 
ted lines) and bulk currents (full lines) computed separately 
with the fitted Rudnick and Stern parameters of panel (c) are 
shown for a harmonic polarization parallel (e) and perpendic- 
ular (f) to the scattering plane. 



the harmonic wavelength (400 nm). Hence, a complete 
theoretical investigation of both intraband and interband 
transitions is required in order to correctly describe the 
nonlinear optical properties of noble- metal nanoparticles. 
In order to go further in the understanding of the bulk 
and surface contributions in the SHG from gold nanopar- 
ticles, the condition d = —b ^ 1 was relaxed, allowing the 
parameter b to take any real values. As shown in Fig. 3c, 
a very good agreement between experimental data and 
FEM simulations is obtained for d — 1, b = 0.1(1) and 
a — 0.5(6)— iO. 2(5). If either the volume contribution (d) 
or the tangential surface contribution (&) is suppressed, 
a less satisfactory fit is achieved as shown in Fig. 3d for 
b = 0: it is not possible to reproduce the plateau ob- 



served at 90 deg nor the intensity maximums reached at 
and 180 deg, indicating that some physics is missing. 
Hence aU contributions are necessary to account for the 
scattered SH intensity, even if in the present case the 
tangential surface contribution is rather weak (see Fig. 
3f). Despite the fact that \a\ < \d\, the normal surface 
current largely dominates the nonlinear response for a de- 
tection polarization perpendicular to the scattering plane 
(Fig. 3e), with an intensity nearly 7.5 larger than the 
other contributions. In contrast, the SHG intensity col- 
lected for harmonic polarization in the scattering plane 
is dominated by the bulk currents as shown in Fig. 3f. 
The ratio between bulk and surface contributions ranges 
from 2 to 6.5 depending on the fundamental polarization 
angle. This clearly underlines the importance of the bulk 
source that can be the dominant contribution in the SHG 
from gold nanostructures depending on the experimental 
configuration!^ 



tion can be discriminated using an interference effect be- 
tween the selected dipolar and octupolar plasmon modes, 
in a specific scattering and polarization configuration. 
By fitting the experimental data by the simulated elec- 
tric fields, we report the first quantitative determina- 
tion of the nonlinear sources efficiencies in spherical gold 
nanoparticles. We show that the actual theoretical mod- 
els, namely the hydrodynamic model and the density 
functional approach, cannot account for the experimen- 
tal data obtained on 150 nm spherical gold nanoparticles. 
This is attributed to the interband transitions that are 
excited resonantly in the present case and not included in 
these models. Finally, we demonstrate that the relative 
local surface and non-local bulk contributions strongly 
depend on the experimental configuration. 



III. CONCLUSION 

In conclusion, it is shown that the local surface and 
non-local bulk contributions to the harmonic genera- 
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